Confinings of QCD at purely classic levels 
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It is well known that confinings and asymptotic freedom are properties of quantum cliromo- 
dynamics (QCD). But tliey can also be observed at purely classic levels. For this purpose we need 
to find solutions to the colorly-sourceful Yang-Mills equations with both confining and asymptotic 
freedom features. We provide such a solution in this paper which at the near-source region is of serial 
form, while at the far-away region is expressed through simple fundamental functions approximately. 
From the solution, we derive out a non-perturbative beta function describing the running of effective 
coupling constant, which is linear in the couplings both in the infrared and ultraviolet region. 
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As is well know, exact solutions to the Einstein equa- 
tion shape our knowledge structure of general relativi- 
ties remarkably. So finding exact solutions to the clas- 
sic Yang-Mills equation like the one of Schwarzschild or 
Riessner-Nordstrom to Einstein equations must also be 
meaningful. But to our surprise, when we search on the 
internet, we almost find no such solutions. The liter- 
ature we obtained most closely related to this topic is 
some on the construction of colorless- monopoles lH-0] 
or colorful black- holes 0-0], or non- gravity but horizon- 
carrying Yang- Mills-Scalar black holes Q. The purpose 
of this paper is to provide solutions to the purely classic 
Yang-Mills theories with both confining and asymptoti- 
cal freedom features. 

Another motivation of this work is reference [oj (see 
also for earlier ideas in a totally different environ- 
ment), in which G. Dvali et al observed that by study- 
ing the classic field stimulated by fixed external sources, 
many features of the underlying quantum theories, such 
as the renormalization group structure, the phenomenon 
of dimensional transmutation, running coupling con- 
stant, asymptotic freedom et al can all be exhibited out at 
the purely classic level. G. Dvali et al take the A(/)^(A < 0) 
theory as examples to illustrate these ideals. Obviously, 
it will be more instructive if for Yang-Mills theories we 
can also see asymptotical freedom and confinement et al 
at purely classic levels. 

Let us begin our stories from the basic action describ- 
ing classic Yang-Mills fields coupled with external quarks 
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Using these notations, the covariant Yang-Mills equation 
can be written as 



1 



'9 



-M^gFn + ^[A^.Fn^ 9..:^ (3) 



where \/—g is the root of metric determinant in some 
specific coordinating system. 

For our purpose, we need to find solutions to the 
above Yang-Mills equation with point like static quark 
sources. The corresponding 4-component quark-current 
vector can be written as 



J^' = (Q"r47r5(f), 0,0,0) 



(4) 



By symmetry reasoning, we know that gauge fields stim- 
ulated by this quark should be spherically symmetric. 
According to reference i and i, the most simple color- 
fields with this symmetry could be written as 

A^edx^" = Atdi + Ardr + Aedd A^d(\) (5) 
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-(C-C^)K{r) (6) 



A^ = --[{C + C^)K{r)sme + D co^e] (7) 

By properly choosing the direction of source quarks' 
charge in the color-space, we can always write the 
matrices C and D involved in the above expressions as 
follows (noting we use i denoting N — i) 

L» = diag.{iV-l,Ar-3,--- ,3-7V,l-7V} (8) 
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substituting this ansatz into the covariant Yang-Mills 
equation, we can easily derive out equations satisfied by 
J(r) and K{r) 



r\j" -2JK^ =g^,,Qr5{r) 
r^K" - K{K'^ - 1 - j2) = 
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where Q is now only a number quantifying the amount 
of color-charges, its direction in the color-space has been 
specified when we choose the form of matrices C and 
D. Although this is a sourceul equation, we will neglect 
the source term when solving differential equations. This 
is very similar to the case of solving Poinsson equations 
V^(^ = eQr5{r) to get Column potentials in electrostat- 
ics or solving Einstein equations R^^y — 2^m^ 

i? = to get 

Schwarzschild metric in general relativity. The differ- 
ence is, in Column potentials and Schwarzschild metric, 
we have definite scheme to relate the solution parameter 
with the source chages (masses is looked as charges of 
gravity). While in the Yang-Mills case, due to the fea- 
ture of asymptotical freedom and confinings, we still have 
no appropriate schemes to relate the solution parameter 
J2, K2 in the following to the source charge Q. 

Now let us give out the first two simple but non-trivial 
solution 



{i) J = 0, K = 0, 
(ii) J = ar + b, K = (b^O) 
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Their non-triviality can be looked out from their corre- 
sponding field strength 
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The first solution describes objects with constant color- 
magnetic field, while the second one describes object 
with both constant color-magnetic field and Column-like 
color-electric field. Obviously, the solution parameter b 
is equal to the source charge Q. If we calculate energies 
stored in this two field configurations (see expression ([20)1 
in the following), we will see that both of them contain 
singularity at r — ^ 0, i.e. E (x dr/r'^. This is a sin- 
gularity similar to that of point-like electric charges in 
classic eletrodynamics. 

The previous two solutions are both non-confining. 
Maybe the more attractive solutions should be the fol- 
lowing two-parameters solution family with both asymp- 



totical freedom and confining features 



{iii)J{r) E J„r«, K{r) ' = ' ^ if„r" (15a) 
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With the corresponding the field strength given by 

r V — 1 IK 
{ttt)For = -A^D, Foe = '^{C + C^) (16) 
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Although converges only at a small region around r — 
0, the serial expressions (|15aP exactly satisfy equations 
(jlOap - ()10b| as long as 
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with J2 and K2 being two parameters characterizing the 
color-electric and color-magnetic charges of the source. 
The analytical expressions (|15bp exactly satisfies ea (|10b|) 
but satisfies (|10aP only up to an order r cos[r ln(r)] oscil- 
lating errors. It can be easily proved that J(r) in this 
expression approaches r In r asymptotically. FIG [T] dis- 
plays one of our numerical solutions explicitly. 

From the field strength expression ea (|16p and the serial 
solution (|15ap . we easily see that both the color-electric 

and color- magnetic field have behavior Fqt, ''^°> 
constant instead of the Column-like increasing r~^. This 
is a feature of asymptotical freedoms. Another way to 
see this is using the fact that J(r) j.2 g^^^j thai 



-> has the meaning of effective color-electric 
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charges. The color-magnetic charge has similar proper- 
ties. The confinement of the solution can be observed from 
the corresponding field-configuration energies, which are 
defined and can be calculated Isl as follows 



E= /d^xT™, T^" 
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FIG. 1: (Color online)Upper part, the left panel is the field 
strength of the usual Columb field, the right panel is that of 
color-electric field with running, asymptotical freedom cou- 
pling constant. Downer part is a typical numerical solution 
to the classic Yang-Mills eas (|10p with the boundary condi- 
tions setting as J{r) J2r^ K{r) 1 + isTar^ The 
dashed line in the figure displays a reference function which 



constant. Changing the value 



is approximately oc ^^^^^^ 
of J2, K2 does not lead to the change of J(r)'s asymptotical 
behavior and K{r)'s oscillating feature. Flipping the sign of 
J2 flips J(r)'s sign correspondingly. 



^ f°°^ (rJ'-J)2 + 2j2x2 + 2r2X'2 + (x2_i)2 
h = / dr (20) 
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In the 
J(r) ^ 



last approximation, we used the fact that 
^ rln(r), K{r) — >■ cos[J{r)]/ J'ir) and in- 
dicated contributions only around r — > 00. Obviously, 
this energy is infrared divergent. Since it is the energy of 
chromo-fields stimulated by colored objects, its diverging 
means that such objects cannot exist as isolating objects, 
and thus confining. 

Now, comes our plaguy question, how to relate the so- 
lution parameter J2, if 2 with the source charge Ql Sim- 
ply using Gauss' law and flux integration does not work 
here. Because, due to the running of effective couplings, 
fluxes on various surface enclosing the charge are not a 



constant. While 5-^,^^ and Q always appear as a whole, 
see eq([TO]). However, this running also suggests us using 
renormalization group method to relate J2, if 2 and Q. 
To use this method [6|, we need first define appropriate 
effective coupling constants(in terms of J(r), K{r) func- 
tions) and prove that their runnings have renormalization 
group structures (RGS). As a try, we define 

(r J' - J\=r, = g, [K^ - l)r=ro = zQ (22) 



r.r~J = g,ir)Q, K^-l = g„,{r)zQ 



(23) 



Obviously, eq([22|) implements the goal of connecting the 
source charge Q to the solution parameter J2 and if 2, 
with z a parameter to be determined latter. While geii") 
and gm{f) defines the effective running color-electric or 
color-magnetic couplings. Reversing eqs p2|) we can write 
J^r^ = ECjnQ"^ K2rl = E4nQ"- Substituting them 
into eq l[23|) we will get 
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If the running of coupling constants ((23l) has RGS, then 
the coefficient c^/"[z, r/ro]'s should satisfy the follow- 
ing structure equation(derived along the routine of Q, 
III. A), 



dx 
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^'"y^'^^.lO], x = ln- (26) 



This is a highly non-trivial constraints on the defini- 
tion of effective couplings But regretfully, we 
find except the lowest two level of coefficients Ci™[z,x], 
C2'"^[z,x], all the other coefficients do not satisfy this 
RGS equation. Maybe, J2, if 2 and Q are not connected 
through so simple a relation like (P^. For example, they 
maybe connected by the following relation 

[{rJ'-jf + 2J^K^ 



a{r = ro) 



2r^if2 + (i^2_i)2] =a{r)Q\27ei) 
1, F[J2rlK2r^]=Q (27b) 



with F{x, y) a properly chosen function different from 
a{ra) so that the running of a{r) has exact RGS. 
Assuming that such a function exist, then from the 
basic equation of motion pH)) and the solution we 
will be able to derive out a classic beta function. 
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where e = 2.71828 • • • , while ProductLog\x] is defined as 
the inverse function oi y = xln[x]. In serial form (9| 



ProductLog [x] = In x — In 
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Eq (l28|) is both a good and a bad result. It's good be- 
cause it gives us a non-perturbative beta function at 
strong coupling. Although classic, it has the same scal- 
ing (3{a) cx a as those \lu\ inspired by supersymmetric 
Novikov-Shifman-Vainshtein-Zakharov(NSVZ) instanton 
analysis and supported by lattice results [HI. It's bad 
because at small couplings limit, its scaling in a does 
not match the quantum expression, /3(a)qunt "~^°> Poo;'^- 
As can be easily proved that, the quantum /3 function 
requires a(r) oc [ln(ro/r)]~^, with vq denoting some ul- 
traviolet cut off. Maybe a natural question is, can we 

r — ^0 1 

set J{r),\K{r)\ — 1 cx [ln(ro/r)]^2 and solve the ba- 
sic equations (fTO|) to get a more correct strong coupling 
beta function? The answer is, we cannot do that way. 
Because, in such cases r^J"_^Q has different signs with 
Jr^o, so the equation J" — 2JK^ — cannot be kept 
consistently at the starting point of integration. It seems 
that as long as asymptotical freedom is required, the clas- 
sic /3 "'^"> 4a is an unavoidable conclusion. 

A possible explanation for this discrepancy is as fol- 
lows. The quantum beta function describes the way the 
couplings change ao — > ai when the theory's UV cut off 
is changed Aq — >■ Ai . It is the result of averaging small 
scale degrees of freedom, i.e. in this case Yang-Mills is 
just taken as an effective description of physics. Obvi- 
ously, when Aq = oo (correspondingly a = 0), we do not 
know how /3(a) is. From the viewpoint of effective theo- 
ries, we also need not to know how it is. While the classic 
beta function ([2H| involves no averaging the small scale 
physics, it takes the theory as an ultimate description of 
physics. It is derived from exact solutions of the theory. 
In other words, Yang-Mills in this case is taken as effec- 
tive theory with Aq = c». It is completely possible that 
/3(a)Ao=oo c>c a while /3(a)Ao=finitc oc a^. In this sense, 
the beta function we defined in this paper is different 
from that defined by G. Dvali et al in [6|. 

As conclusions, we construct three solutions to the 
classic Yang-Mills equation. The first two are exact 
but non-confining. While the third one manifests both 
asymptotical freedom and confinings. It is of serial form 
and exact at the near-source region, but approximate and 
analytical at the far-away region. From this solution we 
derive out non-perturbative beta functions describing the 
running of effective coupling constant, which are linear 
in the couplings both in the infrared and ultraviolet re- 



gion. Our result is sensible to understand confinings of 
QCD at purely classical levels and non-perturbatively. 
It may provide new ingredients for studies of QCD us- 
ing gauge/gravity dualities [l^l- More directly, it will 
be helpful for looking for fully analytical solutions to 
the Yang-Mills equation[l3l. On this point, the value 
of our work is similar to that of Julia and Zee H , , which 
makespreparations for Prasad and Sommerfield's discov- 
eries [8|. 

As a discussion, we must point out that, as a non- 
linear differential equation array, pH)) also allows other 
solutions which are different from what we provide in 
this paper. For example, it can be easily verified that 
Jr^o Kr^Q —7' 1 also leads to reasonable solu- 

tions. More solutions are possible if Kr^o -fr 1- We 
numerically find that all such solutions are confining but 
not necessarily asymptotic free. This point to a conclu- 
sion that, confinings have no direct relation with asymp- 
totical freedom, which is just the viewpoint of [llj- The 
real deep reason of confimcnt is the topological structure 
of the gauge group. 



Appendix 

In numeric practices, we find that no matter how we 
set the r — > boundary condition, as long as the equa- 
tion can be smoothly integrated from r = to 
r ^ oo, we always get confining solutions at r — oo. 
So it maybe very interesting if we can prove that J(r) oc 
r In (r) analytically, since it means that we prove the con- 
finings of classic Yang-Mills theory somehow. Let us try 
the following idea. Consider possibilities. 



(29) 
(30) 
(31) 
(32) 



where p is allowed to take any real values. Of course, enu- 
meration of — 1 ~ J^'s asymptotical behavior cannot 
be completely listed. But through careful examinations 
in the following, we will see that for self-consistencies, 
further listing of other possibilities such as e^^ is of no 
use. Substituting these possibilities into equation (jlOb|) . 
we will get 
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bi) K (xCir^I_i[-r^] + C2r^Ii[-r^] 



> r 2 4 g p 



ci) K 
di) K 



(34) 

, ^^cos[rln(r)] ^ sin[r ln(r)] 



gi'ln(r) 



-rln(r) 

C2 . _ (36) 



v/l + ln(r) Vl + ln(7-) 



where I,y{x) and Ji/(a;) denote the usual Bessel functions. 

In case a) , \i p = 2 then both the two equations (|10ap - 
(llObp can be solved exactly with results J cx Cir^ + 
Car-i, /i: cx DiCOsH+DssinH. SoK^-l~J'^ -r'^ + 
0[r], which is obviously inconsistent with the assumption 
([29|) . So we need only to explore the fact that 



if p > 2, 



if p < 2 
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r-i-OO V ^_r(2-p)/4 (37) 

^ 7-4 e^-p 



eqdMI) = 
ea PUal) 



4_^(2-p)/4 (38) 



-r 1 ^OO ^ 

J > r 4 e^-p 



In either condition, the final conclusions from equations 
and pOap are inconsistent. So, the a) case can be 
excluded. Similarly, the case b) can also be excluded. 
But case c) cannot be excluded by self-consistency. While 
in case d) 



if Ci = 0, 



if C2 = 0, 



7-2 ln2(r). 



(39) 



eq lllOap J ''7°°> ?' In r 
ea pOap =^ J 7 > 



So this is also excluded. The only possibility which 
cannot be excluded is c). While through observing 
equations (f29l) . (|37l) and (f38|) . we know that for self- 
consistency, the asymptotic of — 1 — should be 
something between — and —r'^'^^, but cannot be of 
the r'P form. So, up to corrections of order r2ln[r], 
~~r'^\\i'{r) is almost the unique possibility admitted by 
self-consistency. This finishes our provement of confin- 
ings of QCD at the purely classic levels. 
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